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Abstract. 

We obtained BVIq photometry of IRC— 10443 in 85 different nights dis- 
tributed over two years, and in addition low resolution absolute spectro- pho- 
tometry and high resolution Echelle spectroscopy. Our data show that IRC 
— 10443, which was never studied before in any detail, is a SRa variable, charac- 
terized by A£> = 1.27, AV=1.14 and A/=0.70 mag amplitudes and mean values 
<-B>=13.75, <y>=11.33 and </c>=6.18 mag. Two strong periodicities are 
simultaneously present: a principal one of 85.5 (±0.2) days, and a secondary 
one of 620 (±15) days, both sinusoidal in shape, and with semi-amplitudes 
AV=0.41 and 0.20 mag, respectively. IRC-10443 turns out to be a M7III star, 
with a mean heliocentric radial velocity —28 km/s and reddened by =0.87, 
a third of which of circumstellar origin. The same 0.5 kpc distance is derived 
from application of the appropriate period- luminosity relations to both the prin- 
cipal and the secondary periods. The long secondary period causes a sinusoidal 
variation in color of 0.13 mag semi-amplitude in V — Iq, with IRC— 10443 being 
bluest at maximum and reddest at minimum, and with associated changes in 
effective temperature and radius of 85 K and 6%, respectively. This behav- 
ior of colors argues in favor of a pulsation nature for the still mysterious long 
secondary periods in AGB stars. 
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1. INTRODUCTION 

IRC-10443 (= RAFGL 2209 = NSV 11129 = BD-12°5123) is a bright (#=1.8 
mag) infrared source discovered during the Two Micron Sky Survey (Neugebauer 
and Leighton 1969), that lies in the general direction of the Scutum Star Cloud. 
IRC-10443 was detected by the AFGL survey (Price and Murdock 1983) at 
4.2 /jm, and by IRAS satellite at 12 and 25 /j,m. Its 2MASS magnitudes and 
colors are # s =1.92, J — H=1.35, J — #=1.80. Its spectral type is reported to 
be M6 by Neckel (1958) and Hansen and Blanco (1975), and M6.5 by Nassau 
et al. (1956). IRC— 10443 is present in the NSV catalog of suspected variables 
because J-band observations, obtained at five different epochs (from 21-08-1963 
to 28-06-1965) which are reported in the IRC catalog, seem to trace a variation 
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from magnitude 6.4 to 6.9 (however, the uncertainty of the single measurement is 
similar to the dispersion of the five IRC measurements around their mean). Not 
much more is known about IRC— 10443 and its nature. In this paper we report 
on our BVI photometric monitoring (85 nights distributed over two years) and 
optical spectroscopic observations (low and high resolution) of this object, and 
how they constrain its basic properties. 

2. OBSERVATIONS 

2.1. Photometry 

BVIc CCD photometry of IRC— 10443 was independently obtained with two 
separate telescopes: (a) the 0.30-m Meade RCX-400 f/8 Schmidt-Cassegrain tele- 
scope owned by Associazione Astrofili Valle di Cembra (Trento, Italy). The CCD 
was a SBIG ST-9, 512x512 array, 20 ^im pixels =1.72"/pix, with a field of view of 
13'xl3'. The B filter was from Omega and the Vic filters from Custom Scientific; 
and (6) the 0.50-m f/8 Ritchey-Cretien telescope operated on top of Mt. Zugna by 
Museo Civico di Rovereto (Trento, Italy) and equipped with Optec BVIc filters. 
The CCD was an Apogee Alta U42 2048x2048 array, 13.5 fim pixels = 0.70"/pix, 
with a field of view of 24'x24'. 

The comparison star for B and V bands was TYC 5699-6341-1, for which we 
adopted _B=11.016 and y=10.334 in the standard Johnson UBV system. They 
were obtained from Tycho Bt,Vt data following Besscll (2000) transformations. 
The comparison star for Iq band was TYC 5699-6348-1 for which we adopted 
Jc=6.39, V — Jc=0.55 from the Hipparcos catalog. We had no alternatives for the 
comparison stars. In fact, these two are the only stars within the CCD field of 
view of IRC— 10443 that (i) have reference magnitudes available in literature, (ii) 
are bright enough to be well exposed on the single CCD image without risking to 
saturate IRC— 10443, and (Hi) are photometrically stable to better than 0.02 mag. 
The last point was verified by noting that on all frames in any band we obtained, 
the relative magnitude of the two comparison stars was stable at this level. 

All photometric measurements were corrected for instrumental color equations 
derived nightly by observations of Landolt (1992) standard fields. The good consis- 
tency of the data obtained independently with two different instruments reinforce 
our confidence in the accuracy of the results, in spite of the very red colors of 
IRC— 10443, that are not reached by typical Landolt standard stars. Our photom- 
etry is presented in Table 1. It covers the period from 16-07-2006 to 11-07-2008, 
with observations collected in 85 different nights. The Poissonian component of 
the total error budget is less than 0.01 mag for all the data. The r.m.s. of the 
Landolt standard stars around the color equations they contributed to calibrate 
was on the average 0.019 mag for B, 0.022 for V and 0.031 for Iq bands. 

2. 2. Spectroscopy 

A low resolution, absolutely fluxed spectrum of IRC— 10443 was obtained on 
June 24.97, 2008 UT with the B&C spectrograph of INAF Astronomical Ob- 
servatory of Padova attached to the 1.22m telescope operated in Asiago by the 
Department of Astronomy of the University of Padova. The slit, aligned with the 
parallactic angle, projected onto 2 arcsec on the sky, and the total exposure time 
was 1860 sec. The detector was an ANDOR iDus 440A CCD camera, equipped 



Table 1. Our CCD photometry of IRC-10443. ft and b 
identify the telescopes described in sect. 2.1. 
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Days since Jul 17.0, 2006 (HJD = 2453933.5) 

Fig. 1. Light- and color-curves of of IRC 10443. Open circles: telescope o; filled 
dots: telescope b (see sect. 2.1). 



with a EEV 42-10BU back illuminated chip, 2048x512 pixels of 13.5 //m size. A 
300 ln/mm grating blazed at 5000 A provided a dispersion of 2.26 A/pix and a 
covered range extending from 3250 to 7890 A. 

High resolution spectra of IRC— 10443 were obtained on June 10.04 and July 
22.95 2008 UT with the Echelle spectrograph mounted on the 1.82m telescope 
operated in Asiago by INAF Astronomical Observatory of Padova. The detector 
was a EEV CCD47-10 CCD, 1024x1024 array, 13 /im pixel, covering the interval 
3600—7300 A in 31 orders. A slit width of 200 ^m provided a resolving power 
i? P =26 000. 

3. RESULTS 

3.1. Photometric variability 

The light-curve presented in Figure 1 clearly shows that IRC— 10443 is indeed 
variable. The recorded variability amounts to A£?=1.27 (from 14.41 to 13.14), 
AV=1.14 (from 11.89 to 10.75) and A/ c =0.70 (from 6.45 to 5.75), around the 
mean values <£?>=13.75, <V>=11.33 and </c>=6.18 mag. 

The variability is obviously periodic (see next section) and Figure 1 shows that 
the stars gets hotter (bluest V — Ic) at V maxima, and cooler (reddest V — Ic) at 
V minima. This behavior of the color is a distinctive features of stellar pulsation. 
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Fig. 2. Brightness-color diagram for the 
observations in Figure 1. The line is the 
path followed by a black-body, reddened 
by Eb-v=0-S7, varying in radius and 
temperature at constant luminosity. 



In fact, Figure 2 plots the V — 
Ic color against the V magnitude, 
showing the clear correlation between 
them. If we take a black-body, red- 
den it by Eb-v =0.87, scale its flux so 
to match the average V band bright- 
ness of IRC— 10443, and let it varies 
at constant luminosity, we obtain the 
line in Figure 2, which is a reason- 
ably good fit to the observed points. 
This indicates that the variability 
displayed by IRC— 10443 occurs at 
constant luminosity, in the form of 
expansion + cooling and contraction 
+ warming, as expected in radial pul- 
sations. The corresponding variation 
in spectral type goes from M6.6 III 
to M7.5 III as indicated on the right 
ordinate axis of Figure 2. These cor- 
responding spectral types have been 
obtained by integrating the V and Ic 
(1994) spectra of M III giants reddened by E B - V =0.87 
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bands on the Fluks et al 

(the amount affecting IRC— 10443, see sect. 3.4) 

IRC — 10443 appears as a bona fide SRa variable. SRa variables are semi- 
regular late-type (M, C, S or Me, Ce, Se) giants displaying persistent periodicity 
and usually small (Ay<2.5 mag) light amplitudes. Amplitudes and light-curve 
shapes generally vary and periods are in the range 35-1200 days. Many SRa differ 
from Miras only by showing smaller light amplitudes (Whitelock 1996). 

3.2. Multi periodicities 

Two main periodicities are at the same time present in IRC— 10443: a principal 
and larger amplitude variation modulated by a 85.5 day period, and a secondary 
and smaller one of 620 days. The following expression corresponds to the curve 
fitting the T^-band data in Figure 3 (where t is in HJD - 2450000): 



V(t) = 11.38(±0.02) + 0.41(±0.02) sin 



0.22(±0.02)sin 



i-4042.0(±0.3) 
85.5(±0.2) 

t- 4141 (±10) 



(1) 



620(±15) 

The combination of these two plain sinusoids provides a reasonably close fitting 
to the observed light-curve. Nevertheless, the residuals are larger than the obser- 
vational errors, and an additional weaker component (either periodic or irregular) 
is probably present. Our present data are insufficient to characterize such an ad- 
ditional component, and to resolve it a much longer photometric monitoring is 
required, which we plan to pursue. 

3.3. Spectral classification and radial velocity 

The low resolution spectrum we obtained of IRC —10443 was compared with 
the digital spectral atlas of Fluks et al. (1994), that includes all spectral types from 
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Days since Jul 17.0, 2006 (HJD = 2453933.5) 

Fig. 3. Fitting of the V and V — Jc light-curves of IRC- 10443 with 
sinusoids of 85.5 and 620 days periods (see sects. 3.2, 3.6 and Eq. 1 
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MO to M10 with spectra covering the whole optical range. They are of high flux 
accuracy and of a resolution similar to ours. Literature data suggest a M6/M6.5 
spectral type for IRC— 10443, but our spectrum is quite poorly fitted by the M6III 
reference spectrum from Fluks et al. (1994) library, while the match is perfect 
with a M7III spectrum, as shown in Figure 4. In view of the pulsation activity 
present in IRC— 10443, the difference between our and other spectral classifications 
present in literature can be accounted for by the changes in surface temperature 
that characterize the pulsation activity (see right hand-side ordinates of Figure 2). 

Figure 5 displays a portion centered on Ha of our high resolution Echelle 
spectrum of IRC— 10443 for July 22.95, and by comparison those of bracketing 
spectral types from the atlas of Bagnulo et al. (2003), degraded to the resolution 
of our Echelle spectrum via a Gaussian filter. The spectral progression in Figure 5 
confirms the M7III classification for IRC-10443. 

Mira variables displays emission lines, mainly the higher lines in the Balmer 
series, peaking in intensity at maximum brightness for both O- and C-rich varieties 
(e.g. Panchuk 1978, Yamashita et al. 1977, Mikulasek and Graf 2005), with large 
excursion in intensity along the pulsation cycle. When the first Echelle spectrum 
was exposed on June 10, IRC— 10443 was on the rise and close to maximum bright- 
ness (pulsation phase 0.86), while for the July 22 spectrum it was declining and 
close to minimum brightness (pulsation phase 0.37). Both spectra do not show 
emission in the Ha line, in agreement with the fact that the presence of emission 
lines is far less frequent in SR than in Mira variables. 

The radial velocity of the M7III star is -24.1(±0.8) and -31.4(±0.7) km/s on 
the June 10 and July 22 spectra, respectively. The difference is well accounted 
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for by the pulsation activity. The two observations are separated in time by 
exactly half of the main 85.5 day pulsation period, and their mean value —28 km/s 
can be taken as representative of the systemic velocity of IRC— 10443. At its 
galactic coordinates (1=20°, 6=— 3°) and distance (0.5 kpc, see below), the radial 
velocity expected from galactic disk rotation is +7 km/s (cf also Brand and Blitz 
1993). The 35 km/s difference with the observed radial velocity, suggests that 
IRC— 10443 docs not belong to the young disk galactic population onto which it is 
seen projected and instead it is related to an older population. This is confirmed 
by the high tangential velocity, 86 km/s, derived from the proper motion listed in 
the NOMAD catalog (Zacharias et al. 2004) and the distance estimated in sect. 3.5 
below. 

3.4- Reddening 

The fit with the Flucks et al. (1994) M7III reference spectrum presented in Fig- 
ure 4 constrains the reddening affecting IRC— 10443. The best match is obtained 
with i?B_y=0.87±0.02 for a standard i?y=3.1 reddening law. 

The intrinsic B — V color of M giants does not depend from the spectral type 
and hence the effective temperature, as illustrated by Johnson (1966), Lee (1970) 
and Fitzgerald (1970). Their tabular compilations provide <{B — l / ) >=+1.544 
as the mean intrinsic color of M5 to M8 class III giants. The mean B — V color 
of our observations is < (B — V) >=+2.418, which corresponds to a reddening 
£ B _v=0.87±0.04 affecting IRC 10443. 

These two independent methods converge to the same amount of reddening 
affecting IRC 10443, i?B-v=0.87±0.03, which is adopted in this paper. A fraction 
of this total reddening is probably of circumstellar origin, as supported by the 
detection by Kwok et al. (1997) of emission from circumstellar dust in IRAS low 
resolution spectra. IRC— 10443 lies in the general direction of the Scutum Star 
Cloud (SSC, centered at 1=27°, b=— 3° galactic coordinates). SSC is one of the 
regions of the Milky Way with the highest stellar density, caused by unusually low 
extinction over its area. Reichcn et al. (1990) presented the results of a detailed 
investigation of the extinction over the SSC based on ground-based and balloon- 
born UV survey data. Toward the direction to IRC— 10443 they found that the 
interstellar reddening linearly increases with distance until Eb-v^0.55 is reached 
at 0.5 kpc. Longward, the further rise of the reddening with distance is very slow, 
reaching Eb-v^0.65 at 4 kpc. Only for distances d>6 kpc the reddening increases 
to E B -v~l (Madscn and Reynolds 2005). 

Following the results of Reichen et al. (1990), we therefore conclude that ~l/3 
of the total i?B_y=0.87 reddening affecting IRC— 10443, is of probable circum- 
stellar origin. 

3.5. Distance 

In a seminal paper, Wood et al. (1999) used MACHO observations of late-type 
giant variables in LMC to produce a period-luminosity diagram for them, where 
five separate period-luminosity sequences were identified. Comparing with the 
model prediction of Wood and Sebo (1996), three of these sequences were found 
to coincide with the fundamental and first overtones pulsation modes, while the 
other two seemed to trace the variability induced by ellipsoidal distortion of RGB 
and AGB giants harbored in binary systems. 

Since then, the availability of huge sets of data from large micro-lensing sur- 
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4000 5000 6000 7000 A(A) 



Fig. 4. Absolute spectrophotometry of IRC-10443 for 24.97 June 2008 UT 
(thick line) with superimposed the reference spectrum of a M7III star from the 
atlas of Fluks et al. (1994) reddened by £ B _y=0.87 (thin line). 

veys (MACHO, OGLE, EROS, MOA) contributed to rapidly refine the picture, 
with now up to 14 different period-luminosity sequences identified (e.g. Kiss and 
Bedding 2003, Ita et al. 2004, Soszynski et al. 2005, Soszynski et al. 2007). 

The most recent calibration of the various period-luminosity relations for late- 
type giant variables has been presented by Soszynski et al. (2007). Their relation 
for O-rich semi-regular variables of LMC takes the form K s = — 4.35(logP — 
2.0) + 11.25, where K s band is that of the 2MASS survey. Whitelock et al. (2008) 
have shown that any chemical abundance effect on the A"-band period-luminosity 
relation of Miras must be small. Working with the revised Hipparcos parallaxes of 
van Leeuwen (2007), Whitelock et al. (2008) have derived that period- luminosity 
relation of O-rich Miras in our Galaxy has the same slope and it is 0.1 mag brighter 
than the corresponding one for the LMC. We assume that a similar 0.1 mag shift 
would make the Soszynski et al. (2007) relations applicable to the O-rich SRa 
variables of our Galaxy. Adopting this 0.1 mag shift, a LMC distance modulus 
of (to - M) = 18.39 (van Leeuwen et al. 2007), a LMC reddening of E B - V =0.06 
(Mateo 1998), the extinction relation Aks = 0A42Eb-v for an M-type spectral 
distribution and a standard Rv—3.l extinction law (Fiorucci and Munari 2003), 
the distance to IRC— 10443 corresponding to the 85.5 day period is 0.5 kpc. 
Soszynski et al. (2007) relation for the long secondary periods of O-rich red giants 
in LMC takes the form K s = -4.41 (log P - 2.0) + 15.05, and when applied (with 
the same 0.1 mag shift as above) to the 620 day secondary periodicity displayed 
by IRC— 10443, it provides the same distance, 0.5 kpc, as obtained for the 85.5 
day main period. Such 0.5 kpc distance is adopted for IRC— 10443 in this paper. 

3.6. On the nature of the long secondary period 

In spite of large investigation efforts, both observational and theoretical, 11 the 
cause of the long secondary periods seen in cool giants remains a mystery at the 
present time" as recently remarked by Wood (2007). 

It has been known for a long time that some semi-regular variables shows 
the presence of a long secondary period (LSP) in their light-curves, typically ten 
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M6III 



M5III 



times longer than the primary pulsation period. This ratio for IRC— 10443 is 7.25. 
Lists of local giants displaying LSPs have been published, among others, by Houk 
(1963), Mattei et al. (1997), Kiss et al. (1999). Wood et al. (1999) found that 
~25% of all variable AGB star in LMC show LSPs. A similar fraction, ~30%, of 
local semi-regular variables has been found by Percy et al. (2004) to display LSPs. 

The semi-regular variables 
appear to pulsate in the 
first overtone (Lebzelter and 
Wood 2006), and thus it 
could be tempting to relate 
the LSPs with pulsation in 
the fundamental mode. How- 
ever, as found in theoreti- 
cal models by Fox and Wood 
(1982) and verified by ob- 
servations (e.g. Kiss et al. 
1999, Mattei et al. 1997), 
the ratio of fundamental to 
fist overtone periods is close 
to two, ruling out that LSPs 
are due to pulsations in the 
fundamental mode. Estab- 
lished observational facts are 
that LSPs are accompanied 
by radial velocity variations 
(of a lower amplitude than 
observed for the primary pe- 
riod; Hinkle et al. 2002, 
Wood et al. 2004) and by 
variation in intensity of the 
Ha absorption (that could 
trace a variable filling by an 
emission component of chro- 
mospheric origin; Wood et al. 
2004). In addition, cool variable giants showing LSPs rotate at similar velocity 
and show similar dust-free IRAS colors as cool variable giants not showing the 
LSPs (Olivier and Wood 2003). In RGB objects showing LSPs, the period of asso- 
ciated radial velocity variations is twice the period of photometric LSP variability 
(as expected in the case LSPs arise from ellipsoidally distorted giants in binary 
systems; Adams et al. 2006), while in AGB objects the period is the same (as 
expected in the case of pulsations; Wood et al. 2004). Various explanations of the 
LSP phenomenon have been proposed, but all have encountered some problems, 
as discussed by Wood (2007, and references therein). 

A striking feature displayed by IRC— 10443 is illustrated in the bottom panel of 
Figure 3, where the color variation is fitted with two sinusoids of the same periods 
of those fitting the V light-curve (cf Eq. 1 and the top-panel of Figure 3). Their 
semi-amplitudes are 0.23 mag for the sinusoid associated to the principal 85.5 day 
period, and 0.13 mag for the LSP, with a mean V — Jc=5.19. From Figure 3, it 
is evident that IRC— 10443 is bluest when it is at the maximum brightness along 



6550 



6560 



6570 6580 



Fig. 5. Portion, centered on Ha and Cal 6572.8, 
of the high resolution 22.95 July 2008 Echelle 
spectrum of IRC — 10443. Spectra of cool giants 

from Bagnulo et al. (2003) are plotted for 
reference. All spectra are continuum normalized, 
offset in ordinates for better visibility, and shifted 
to 0.0 radial velocity. 
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the LSP cycle, and reddest when it is at the minimum brightness. This is the 
same pattern observed for the principal 85.5 day period and strongly argues in 
favor of a pulsation interpretation of LSP phenomenon, at least in IRC— 10443. 
The simultaneous presence of two pulsations also accounts for the dispersion of 
the points in Figure 2 along the back-body curve. The dispersion would have been 
significantly reduced if only one pulsation would have been present, as confirmed 
by removing from observations one or the other of the sinusoidal variation of the 
color and re-plotting Figure 2. 

The variation in V — Iq color can be transformed into variation in effective 
temperature of the underlying stellar photosphere using the reference continuum 
energy distribution given by Fluks et al. (1994) along the spectral sequence of 
M giants. The total amplitude of 0.46 mag observed in V — Iq for the principal 
period, corresponds to a change in 0.73 spectral types around the M7III mean, 
and therefore to a variation from 3030 to 3175 K in effective temperature. The 
0.26 mag total color amplitude of the LSP would correspond to a change of 0.41 
spectral types around the M7III mean, meaning a variation from 3060 to 3145 K. 
If both pulsations are supposed to occur at constant luminosity, the corresponding 
total excursion in radius is ~10% for the principal 85.5 day period, and ^6% for 
the LSP. 

These excursions in radius are less than that inferred by Wood et al. (2004) 
from radial velocity observations at optical wavelengths of a sample of three hot- 
ter AGB stars characterized by longer LSP than IRC— 10443. We estimated the 
change in effective temperature and radius of the underlying photosphere. If in- 
stead we had referred to a black-body fitting of the observed optical spectrum 
of the star (re-shaped by the extremely strong TiO molecular absorptions), the 
variation in color temperature and radius of IRC— 10443 would have been almost 
twice larger, 150 K and 10 %, respectively. 

It is worth noticing that recent theoretical improvements in the treatment of 
pulsation, like inclusion of time dependent turbulent convection (Olivier and Wood 
2006), are opening new modeling possibilities for pulsation modes in cool giants. 
Important applications to the long lasting problem of what is driving the LSPs 
could be obtained in the near future (cf. Wood 2006). To better characterize the 
object and increase its interest as a test target for current theories, we plan to 
continue a tight observational monitoring of IRC— 10443 over the next years, long 
enough to cover at least the whole next LSP period. 
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